*Import the newly registered data for the years 2016-2019: 

clear all

local count=1 

forvalues num=2016/2018 {

frame create cars_`count'
frame change cars_`count'

import delimited "$root/Data/Original/NEUZU-`num'.txt"

keep if fahrzeugart_code==001

rename erstinvekehrsetzung_kanton erstinverkehrsetzung_kanton

keep typengenehmigungsnr treibstoff_code erstinverkehrsetzung_kanton erstinverkehrsetzung_monat plz erstinverkehrsetzung_jahr marke_und_typ 

gen month_year=ym(erstinverkehrsetzung_jahr, erstinverkehrsetzung_monat)
format month_year %tm
gen quart_year=qofd(dofm(month_year))
format quart_year %tq

local count=`count'+1
}

forvalues num=2019/2022 {

frame create cars_`count'
frame change cars_`count'

import delimited "$root/Data/Original/NEUZU-`num'.txt"

keep if fahrzeugart_code==001

keep fahrzeugart_code typengenehmigungsnr treibstoff_code erstinverkehrsetzung_kanton erstinverkehrsetzung_monat plz erstinverkehrsetzung_jahr marke_und_typ 

gen month_year=ym(erstinverkehrsetzung_jahr, erstinverkehrsetzung_monat)
format month_year %tm
gen quart_year=qofd(dofm(month_year))
format quart_year %tq

local count=`count'+1

}

frame change cars_1 

forvalues i=2/7 {
	
	frameappend cars_`i', drop 
}


gen BE=(erstinverkehrssetzung_kanton=="BE")

preserve

gcollapse (count) freq_reg=typengenehmigungsnr (sum) reg_BE=BE, by(typengenehmigungsnr year_reg)

rename typengenehmigungsnr tgcode

*store the data: 
save "$root/Data/Original/NEUZU_CH.dta", replace

restore 

**Create the NEUZU_TOT database that is used for the EV subsidy description in the appendix: 
frame put Fahrzeugart_Code tm reg_canton, into(ovr_`num')

rename erstinvekehrsetzung_kanton reg_canton

gcollapse (count) nb_cars_= fahrzeugart_Code (firstnm) quart_year, by(treibstoff_code month_year reg_canton)

reshape wide nb_cars , i(month_year reg_canton) j(treibstoff_code) string

egen nb_cars_tot=rowtotal(nb_cars_*)

save "$root/Data/Original/NEUZU_TOT.dta", replace
